Propagation and Ghosts in the Classical Kagome Antiferromagnet 
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We investigate the classical spin dynamics of the kagome antiferromagnet by combining Monte 
Carlo and spin dynamics simulations. We show that this model has two distinct low temperature 
dynamical regimes, both sustaining propagative modes. The expected gauge invariance type of the 
low energy low temperature out of plane excitations is also evidenced in the non linear regime. A 
detailed analysis of the excitations allows to identify ghosts in the dynamical structure factor, i.e 
propagating excitations with a strongly reduced spectral weight. We argue that these dynamical 
extinction rules are of geometrical origin. 

PACS numbers: 75.10Hk,75.40Gb,75.40Mg,75.50.Ee 



Geometrical frustrated magnets are currently a source 
of high interest for the exotic phases and unexpected dy- 
namics that they are liable to generate. A full insight 
about their behaviors is still far from having been ac- 
quired, in particular at the lowest temperatures. 

A prototype is the classical Heisenberg kagome anti- 
ferromagnet As a basic distinctive feature of the ge- 
ometrical frustration, its ground state consists in a con- 
tinuous connected manifold of spin configurations. At 
high temperatures (T/ J > 0.1, with J the first neighbor 
exchange), the system is paramagnetic. It enters what 
we shall call from now on a cooperative magnetic phase 
in the range 5 • 10"'^ ^ T/J ^ 0.1 where short range 
correlations are enhanced. At the lowest temperatures 
(T/ J < 5 • 10~^), thermal fluctuations above each of the 
spin configurations of the ground state manifold are not 
equivalent and drive an entropic based order out of dis- 
order mechanism ^ , eventually selecting a spin plane Q 
and developing an octupolar order Q . We shall call this 
phase coplanar to distinguish it from the former. While 
in both low temperature regimes it was shown that spin 
pair correlations remain short ranged it is only in 
the coplanar phase that the continuous degeneracy of the 
manifold was argued to be reduced to a discrete one, de- 
scribed by the 3-colorings of the lattice 0]. Altogether, 
these results provide a rather clear picture of the thermo- 
dynamics of the classical kagome antiferromagnet, which 
should apply to experimental compounds with large mag- 
netic moments but also be of some relevance for low spin 
systems, since quantum fluctuations often play a signifi- 
cant role at very low temperatures only. 

A much poorer understanding of the spin dynamics is 
in contrast available. To our knowledge, only one nu- 
merical investigation was so far conducted 0] , which fur- 
thermore was not resolved in momentum vectors Q, thus 
ignoring any diffusive or propagating aspects of the ex- 
citations. In this letter, we analyze the temperature de- 
pendent dynamics of the classical kagome antiferromag- 
net from two point of views. We first show that at low 
temperatures, spin waves (SW) do propagate and are sen- 



sitive to the underlying spin texture, either cooperative 
paramagnetic or coplanar. Quantitative analysis of the 
dynamical structure factor is performed and provides the 
characteristic time scales. Additionally, the invariance of 
the linear SW spectra with respect to the ground state 
spin configurations on which they are built is evidenced 
in a wide range of temperatures, including those where 
non linear effects are at play. We next put forward that 
peculiar excitations develop that would be almost invis- 
ible to dynamical spin-pair correlations sensitive probes, 
such as inelastic neutron scattering. 

The numerical method used in this work is a combi- 
nation of an hybrid Monte Carlo (MC) method, which 
allows generating samples of spin arrays at a given tem- 
perature, and an integration of the non-linear coupled 
equations of motion for the spin dynamics (SD): 
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where j is a first neighbor of i and J > is the antifer- 
romagnetic exchange Q. The numerical integration has 
been performed up to < = 1000 J^^ using an 8'^-order 
Runge-Kutta method (RK) with an adaptative step-size 
control. The RK error parameter as well as the RK or- 
der have been fixed in order to preserve the euclidian 
distance with a test-full diffusion of Eq. [1] performed 
with the more robust but time consuming Burlisch-Stoer 
algorithm. As a result, trivial constants of motion, such 
as the total energy Etot and magnetization Mtot, are con- 
served with a relative error smaller than 10~^. As for the 
spin arrays samplings by the MC method, a first run has 
been performed in order to find an optimal set of tem- 
peratures for a parallel tempering scheme Q , which min- 
imizes the ergodic time [l3] ■ A reduction of the solid an- 
gle for each spin fiip trial together with rotations around 
the local molecular fields ensure a rate acceptance above 
40%. The numerical simulations reported in this work 
were performed on samples of L x L x 3 spins with L = 36 
and periodic boundary conditions. Our interest lies in the 
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FIG. 1: (color online) Intensity map (a. it.) of the scattering 
function vs. a; and Q = [h, 0), for T/J = 5 • lO-"*. The LSW 
dispersion relations t<;(Q) for the phase q = \/3 x are 
plotted as blue lines. The white dotted line corresponds to 
the constant-Q scans presented on Fig. [21 



scattering function, namely the time and space Fourier 
transform of the dynamical spin-pair correlations: 



dt 



2t:N 



(S,(0)-S,(t))e 



— 'iQ-R,ij —iLut 



(2) 



where Q and uj are the momentum vector and energy 
transfer, (...) is the ensemble average, 
and N the number of spins. 

Details about static properties will not be given here. 
It however is worth noting that our results for uj ~ 
at very low temperatures {T/J < 5 • lO"-^), when en- 
tropic selection is at work, points towards fluctuations 
predominantly associated with the so-called q ~ a/3 x a/3 
phase. This meets with previous conclusions, although 
those were inferred from instantaneous ensemble aver- 
ages [1, @, [lH. We similarly got very good agreements 
with previous numerical or analytical investigations of 
the specific heat, the coplanar ordering, or the instanta- 
neous scattering function [1, 0, 0, , giving confidence 
on the quality of our numerical simulations. 

Let us now focus on the dynamical properties of the 
kagome antiferromagnet. Although the propagation of 
collective excitations may appear unexpected in such a 
system, where the spin-pair correlation function decays 
exponentially with distance at finite temperatures 0], a 
sufficient temporal and spatial stiffness may lead to the 
propagation of SW in locally ordered regions. There- 
fore, a required condition for the development of SW 
excitations is an increase with decreasing temperature 
of the autocorrelation time Ta featuring the lifetime of 
locally ordered states. Tq has been numerically evalu- 
ated by integrating the scattering function over all Q- 
values in the reciprocal space, which gives access to the 
time Fourier transform of the autocorrelation function 
A{t) = (Si(0) • Si{t)). A fit of the obtained quasiclas- 



tic (QE) signal using a lorentzian shape jif^^, associ- 
ated with a decaying exponential law A{t) — exp (—Fat) 
in time space, allows extracting the Half Width at Half 
Maximum (HWHM) Fq cx 1/T(j. As the temperature is 
decreased, this shows an algebraic variation F^ = AT"^ 
with C = 0.995 ± 0.018 for T/J < 0.1 (see Fig. [SJi), 
transposing to a slowing down of the spin fluctuations, 
but nevertheless no spin freezing even at temperatures as 
low as T/J = 5 • 10~^. Interestingly, the same thermal 
variation is observed in the cooperative paramagnetic 
(T/J > 5 • 10-3) and the coplanar states {T/J < 5 • lO'^) 
regimes, asserting that the entropic selection favoring the 
coplanar manifold has no influence on the lifetime Tq of 
locally ordered states. Now that we have characterized 
the temporal stiffness associated with the T~^ slowing 
down of Ta, we address the question of well defined ex- 
citations as well as their possible propagation. An ev- 
idence of the existence of SW-type excitations at low 
temperatures is explicit in the excitation spectrum for 
T/J = 5 • 10"^ (see Fig. [T]). For comparison, the linear 
spin wave (LSW) spectrum [ll| emerging from the pure 
q = \/3 X a/S phase is shown. The SD simulations evi- 
dence a large weight of S'(Q,w) at this LSW spectrum, 
confirming that the q ~ \/3 x ^/S short range dynamical 
correlations are favored at very low temperature. 

The analysis of the spectrum as a function of the tem- 
perature allows getting more insights about the forma- 
tion of the SW excitations. We show in Fig. [5] the 
frequency dependence of S{Q,uj) at the point Qq = 
27r(3/4, 0) in the reciprocal space, located between the 
Brillouin zone (BZ) boundary and the BZ center (sec Fig. 
[1]), where the soft, acoustic and optical modes arc par- 
ticularly easy to distinguish. These constant-Q scans are 
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FIG. 2: (color online) Temperature weighted scattering func- 
tion vs. energy at different temperatures for Qo — 27r(3/4, 0). 
Inset : position of soft and acoustic modes versus tempera- 
ture; error bars obtained from several fitting processes. 
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represented for temperatures from T/ J = 0.5 to 5 • 10~^. 
At high temperatures (T/J > 0.2), only a QE signal 
centered at a; = contributes to the scattering func- 
tion S{Cl,(jj). A single broad excitation at finite energy 
comes into sight on decreasing T/J from 0.2 to 10~^, al- 
though strongly softened compared to the LSW theory 
expectation (see inset of Fig. [2]). Below T/J = 10~^, 
the broad peak splits into two excitations, respectively 
associated with the SW acoustic modes and the emerg- 
ing soft modes, gradually separating from each other and 
getting thinner as the temperature goes down (see Fig. 
[2]). The softening of the modes dies away to disappear 
below T/J = 5 • 10~^. The soft mode, expected to be 
non-dispersive in LSW theory, is here observed at finite 
energy, due to the non linear nature of Eq. [U which 
takes account of the interactions between the SW. This 
effect is expected to decrease with temperature, which is 
consistent with the fact that the soft mode drops to zero 
energy when temperature goes down (inset of Fig. [2]). 
Finally, one can discern an additional peak at a; ~ 2 J for 
T/J < 2 • 10^^, corresponding to optical modes. 

Each mode i of the excitation spectrum can be char- 
acterized by its dispersion relation a;*(Q), its lifetime 
'^sw (-l^sw)""^' its intensity Iq, all these quantities 
being accessible by fitting the excitation spectrum at dif- 
ferent temperatures and Qo values. Assuming Lorentzian 
shape for magnetic excitations, the scattering function 
writes : 

^^"■"> = I^ (r.p + (r..(Q„)P 

where i runs over soft, acoustic and optical magnetic 
peaks for a particular Qo value. We show Fig. [HJd the 
thermal variation of the resulting SW HWHM Fsw oc 
for Qo = 27r(3/4, 0). It is found out, contrarily to 
Ta (see Fig. [3^), that tsw follows two distinct regimes 
below and above T/J ^ 5 ■ 10~^, both consistent with 
an algebraic law rsw = AT~^ . rsw is in principle re- 
duced by two physical processes. The first, common to 
all magnetic systems, is associated to the thermal fluctu- 
ations and anharmonic interactions between SW modes. 
In a disordered medium, this process is overwhelmed by 
a second one, induced by the motion of the system be- 
tween the different ground states. In other words, even in 
the linear approximation, one is left with a set of linear 
equations of motions with time dependent initial condi- 
tions, the time variations of those being set by the au- 
tocorrelation time Ta- In the pyrochlore antiferromag- 
net, it has been shown that Ta is also proportional to 
T~^ and based on the above interpretation, it was pro- 
posed that the SW lifetime rsw is proportional to T~^/^ 
[l3|. This behavior is expected at least in the cooper- 
ative paramagnetic regime. In the corresponding tem- 
perature range, we find C = 0.18 ± 0.07, which is much 
lower than 1/2. In the low T regime, due to the se- 
lection of coplanar spin configurations, the out-of-plane 
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FIG. 3: (color online) (a) Ta oc obtained by fitting the 
QE signal S{ljj), shown in the inset at several temperatures 
from T/J = 0.5 to 5 • 10"''. (b) Tsw oc t^-^ for the soft 
(blue squares), in-plane (black diamonds) and out-of-plane 
(red circles) acoustic modes for Qo = 27r(3/4, 0). The error 
bars have been obtained by averaging over difi'erent fitting 
processes. If not seen, they are smaller than the symbols. The 
inset displays Tsw as a function ujo — uj{Q), le for several Q- 
values. The hatched region forbids propagating excitations. 



a;^(Q) and the in-plane ti^"(Q) modes become different, 
which allows distinguishing the corresponding scattering 
process S'(Q,cj) = S^iQ,Lu) + S'll(Q,cj). Within the 
LSW, or equivalently at very low temperatures, the out 
of plane scattering function 5'-'-(Q,w) is gauge invariant 
like, i.e does not depend on the three coloring state on 
top of which the excitations develop (see Fig. |4)3, left). 
Conversely, the in-plane contribution S'" (Q, uj) differs for 
each configuration it is build on (see Fig. [4|d, right). 
In this coplanar regime, rsw seems to behave in similar 
ways for the in-plane {( = 0.69 ± 0.22) and out-of-plane 
(C = 0.783 ± 0.057) acoustic modes, whereas ( is slightly 
weaker for the soft mode (C = 0.577 ±0.09). Thereby, we 
clearly see two distinct dynamical regimes with a C value 
in the coplanar phase significantly larger than the value 
in the cooperative paramagnetic phase. This suggests 
that the SW lifetime is sensitive to the entropic selec- 
tion of the coplanar manifold, that latter inducing some 
kind of stiffness in the spin texture. We can also wonder 
about the propagation of these magnetic excitations. If 
Fsw/wo < 1 for a given mode, it can be considered as 
a propagative SW, its lifetime rsw being longer than its 
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FIG. 4: (color online) (a) Intensity map (a. it.) in reciprocal 
space for ui — J and T/J = 5 • 10~*. The first and neighbor- 
ing BZ are in blue, (b) Out-of-plane (_L) (left) and in-plane 
(II) (right) components for q = x \/3 and g = spin 
configurations, (c) Anisotropy parameter R = Imax/Imin vs. 
temperature. Inset : scattering function plotted for different 
temperatures in function of Lp, as defined in Fig. (a). 



period lUq ^. This actually is the case for the soft and 
acoustic modes, at least up to T/J < 0.1 (see inset of 
Fig. 3b for T/J = 5 • IQ-^). This underlines that SW 
propagation in the kagome antiferromagnet, although of 
different nature, is possible in both regimes of coopera- 
tive paramagnetism and entropy induced coplanarity. 

We finally focus on the spectral weight distribution 
in reciprocal space, which is non-uniform for the excita- 
tions emerging from the BZ centers, i.e the out-of-plane 
acoustic like modes. A two dimensional intensity map 
in reciprocal space is shown in Fig. 0^, for uj = J and 
T/J = 5 • 10~^. The spectral weight S (Q, J) reaches its 
maximum value /max along a*, b* or a* — b* axis, and 
fades out in each corresponding perpendicular direction 
(with an intensity /min). This results in the presence 
of "ghosts" in the excitation rings, i.e existing excita- 
tions with a strongly reduced cross section, that would 
be invisible e.g in neutron scattering experiments. Pa- 
rameterizing these rings by the angle tp and integrating 
over a small width 6 of the ring (see Fig. UK) allows to 
quantitatively analyze the spectral weight anisotropy as 
a function of the temperature. Fig. Hj; displays the evo- 



lution of the anisotropy parameter R = Imax/ Imin with 
the temperature. Strong discrepancies between the co- 
operative paramagnetic regime (i? ~ 1) and the entropy 
driven coplanar regime, in which the anisotropy strongly 
increases with decreasing temperature, are evidenced. At 
the lowest temperatures, where SW propagates onto a 
disordered manifold, the strongly fluctuating spin tex- 
ture could therefore be expected to drive these extinc- 
tions. Actually, two distinct arguments rather support 
a purely geometrical origin. Firstly, out of plane excita- 
tions arc gauge invariant like. Therefore, the fluctuating 
nature of the manifold should not play any role. Sec- 
ondly, we have numerically computed (Q, w) for config- 
urations prepared in slightly distorted ordered q = and 
q = \/3 X \/3 phases and performed a LSW expansion 
aroimd these two phases (Fig. [3]d). All calculations re- 
produce this spectral anisotropy, pointing out that the 
ghost excitations rather originate from the peculiar ge- 
ometry of the lattice revealed by spin coplanarity. 

In conclusion, the propagation of spatially structured 
collective excitations has been numerically evidenced and 
quantitatively studied in the classical kagome antiferro- 
magnet. Although the SW exist in both cooperative 
paramagnetic and coplanar regimes, their lifetime was 
found very sensitive to the entropic selection occurring 
below T/J = 5 • 10""^, in contrast with the same inverse 
temperature dependence of the autocorrelation time in 
both regimes. At very low temperatures, these propaga- 
tive modes possess a noteworthy non uniform spectral 
weight expressing effective dynamical extinction rules. 
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